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I.  A *ST* ACT 


Intensity  interferometry,  as  developed  by  Hanbury  Brown  and  Twiss  for  stellar  observ¬ 
ation,  has  shown  relative  insensitivity  to  atmospheric  scintillation.  However,  with 
classical  sources,  the  limitations  placed  on  this  technique  by  quantum  noise  and 
detector  efficiency  are  somewhat  severe.  This  situation  is  vastly  improved  when  larer 
illumination  is  employed.  Generalizing  a  methodology  of  Marchand  and  Wolf ,  the  far- 
zone  behavior  of  the  mutual  intensity  function  is  derived  for  an  intermediate  time 
average.  This  result  is  used  to  reconstruct  the  irradiance  distribution  of  a  spatiall; 
incoherent  source.  The  far-r field  intensity  distribution  is  recorded  spatially  for  ore| 
time-resolution  unit  of  the  detector.  The  resulting  spatial  signal  is  cress  corre-  jj 
lated  with  itself  and  related  to  the  intensity  distribution  over  the  source.  Thus  <j 
without  averaging  in  the  time  domain,  a  spatial  Fourier  transform  relation  is  derived  I 
between  the  far- field  intensity  correlation  and  the  source  irradiance,  similar  to  the  I 
results  of  Hanburv  Brown  and  Twiss.  1 
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INTENSITY  INTERFEROMETRY  IN  THE  SPATIAL  DOMAIN 


ABSTRACT 

intensity  interferometry,  as  developed  by  H anbury  Brown  and  Tviss 
for  stellar  observation,  has  shew  relative  .insensitivity  to  atmospheric 
scintillation.  However,  with  classical  sources,  the  limitations  placed 
on  this  technique  by  quantum  noise  and  detector  efficiency  are  somewhat 
severe.  This  situation  is  vastly  improved  when  laser  illumination  is 
employed.  Generalizing  a  methodology  of  March aud  and  Wolf,  the  far- zone 
behavior  of  the  mutual  intensity  function  is  derived  for  an  intermediate 
tint  average.  This  result  ie  used  to  re  cent  tract  the  irradiance  distri¬ 
bution  of  a  spatial  ly-in  coherent  source.  Hie  far- fie Id  intensity  dis¬ 
tribution  is  recorded  spatially  for  on«  time-resolution  unit  of  the  de¬ 
tector.  The  resulting  spatial  signal  is  cross  correlated  with  itself 
and  related  to  the  intensity  distribution  over  the  source.  Thus  without 
-leveraging  in  the  time  domain,  s  spatial  Fourier  transform  relation  is 
derived  between  the  far- fie Id  intensity  correlation  and  the  source 
irradiance,  similar  to  the  -results  of  Hanbury  Brown  and  Twiss. 
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I.  INTRODUCTION 


The  history  of  intensity  interferometry  is  rooted  in  the  work  of 

i* 

R.  Hanbury  Brown  and  R.  Q.  'tyiss.  Their  earliest  investigations  dealt 
with  the  problem  of  resolving  stellar  radio  sources  by  a  technique  in¬ 
volving  the  correlation  of  the  squared  outputs  of  two  receivers.  The 
advantages  accrued  relate  to  the  reduction  of  certain  kinds  of  experi¬ 
mental  constraints  as  well  as  the  comparative  insensitivity  of  the 

2 

method  to  almospneric  scintillation.  A  preliminary  conclusion  readied 

at  that  time  was  that  this  technique  of  intensity  correlation  would  not 

be  applicable  at  optical  frequencies  due  to  limitations  imposed  by 

3 

photon  noise.  However,  in  later  work,  Hanbury  Brown  and  Twiss  showed 
that  meaningful  intensity  correlations  could  be  made  at  optical  frequen¬ 
cies  even  with  highly  degenerate  sources.  The  limitations  placed  on  this 
approach  by  quantum  noise  and  detector  efficiency  have  been  somewhat 
severe,  calling  for  highly  refined  experimental  technique. 

For  laser  illumination,  the  situation  is  very  different.  The  signal- 

4 

to-noise  ratio  can  be  typically  increased  by  six  orders  of  magnitude. 
However,  the  statistics  of  the  source  must  be  considered  in  the  measure¬ 
ment.  The  key  to  relating  intensity  correlations  to  some  property  in¬ 
volving  field  correlations  lies  in  the  assumption  of  gaussian  statistics,5 
for  which  all  higher  moments  are  determined  from  the  first  and  second. 
Single  mode  lasers,  though,  are  distinctly  non-gaussian  in  their  statis¬ 
tics,  and,  therefore,  cannot  be  described  by  theory  framed  for  thermal 
sources.  Put  with  the  addition  of  only  a  few  axial  modes,  the  field 
amplitude  becomes  nearly  gaussian  distributed.6 

7 

The  prin  up-il  formula  used  by  Hanbury  Brown  and  TViss  to  infer  the 
diameter  of  a  listant  source  shows  that  the  time-averaged  correlation  of 
intensities  at  two  points  is  equal  to  the  product  of  a  function  involving 


* 

References  may  be  found  on  page  29. 
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the  temporal  characteristics  of  the  source  with  the  square  of  the  spatial 
Fourier  transform  of  the  source  intensity  distribution.  If  the  intent  of 
an  intensity  correlation  experiment  is  to  gain  information  concerning  the 
source  intensity  distribution,  then  the  source  temporal  statistics  may  be 
of  little  interest  in  themselves.  Their  consideration  is  necessary  if 
the  product-output  of  the  detectors  is  averaged  in  the  time  domain  (as 
it  nearly  always  is)  to  overcome  the  limitations  imposed  by  photon  and 
detector  noise  and  possibly  reduce  scintillation  effects  produced  by 
transmission  in  the  atmosphere. 

Intensity  interferometry  can  be  understood  as  a  two-point  corre¬ 
lation  of  intensities  following  the  squaring  of  the  electric  field  at 
the  detector.  If  the  source  is  quasi-monochromatic,  each  differential 
element  on  the  object  emits  a  number  of  temporal  inodes  which  interfere 
with  each  other  at  the  detector.  Assuming  the  detector  has  sufficient 

g 

speed  to  detect  these  beat  frequencies,  it  is  the  amplitude  and  phase 
of  the  incoming  intensity  envelope  which  is  utilized.  It  is  then  commonly 
argued  that  the  random  fluctuations  in  the  temporal  statistics  from  dif¬ 
ferent.  points  of  the  source  cause  the  beat  frequencies  from  each  source 
point  to  add  on-the-average  incoherently  at  the  detector.  The  time- 
averaged  intensity  correlation  is  then  proportional  to  the  squared 

spatial  Fourier  transform  of  the  source  intensity.  This  approach  results 

9 

in  essentially  a  squared  version  of  the  Van  Cittert-Zernike  theorem. 

We  wish  to  argue  that  the  requirement  of  surface  roughness  at  the 
source  (to  assure  spatial  incoherence)  is  sufficient  to  guarantee  the 
incoherent  addition  of  beat  frequencies  at  the  intensity  detectors. 

Thus  if  temporal  noise  (photon  noise,  time-dependent  detector  noise)  is 
largely  absent  in  a  local  spatial  sense,  as  might  be  the  case  with  a 
multi-axial  mode  laser  used  with  photographic  detection,  then  the  inten¬ 
sity  information  might  be  gathered  during  one  resolution  time  of  the 
detector  over  a  plane  section  normal  to  the  direction  of  light  propa¬ 
gation.  Any  noise  arising  in  the  process  would  be  spatial  in  nature, 
and  might  be  averaged  out  by  taking  a  sufficiently  large  area  of  spatial 
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correlation.  The  reduction  of  atmospheric  spatial  noise  would  be  similar 

10 

to  a  process  known  as  aperture  averaging.  Film-grain  noise  would  be 
extremely  well  averaged  due  to  the  relatively  large  area  of  averaging. 

The  relative  insensitivity  of  intensity  interferometry  to  turbulence 
is  due  to  the  assumed  dispersionless  nature  of  the  propagation  medium. 
Since  each  temporal  frequency  sees  the  same  refractive  index,  the  dif¬ 
ferential  (beat)  frequencies  remain  unchanged.  However,  the  spatial 
Fourier  transform  relation  between  the  source  and  the  far- field  scales 
as  the  average  frequency,  not  the  beat  frequency,  end  thus  the  resolution 
afforded  by  optical  frequencies  is  maintained. 

The  concept  of  examining  spatial  beat  frequencies  of  second-order 
correlation  is,  of  course,  not  new.  Many  classical  fie id- correlation 
interferometers,  as  well  as  holographic  experiments,  are  built  on  this 
principle  involving  a  spatial  or  time  lag  between  interfering  beams  of 
the  same  source.  Mare  difficult  is  the  spatial  recording  of  beats  from 
two  independent  sources  as  demonstrated  by  Magyar  and  Mandel.*1 

We  have  examined  the  problem  of  relating  far- field  spatial  intensity 

(fourth-order)  correlations  to  the  intensity  distribution  of  the  source, 

without  here  considering  the  limitations  due  to  noise.  Our  approach  has 

Deen  completely  classical,  drawing  on  a  straightforward  generalization 

12 

of  a  methodology  given  recently  by  Marchand  and  Wolf.  Our  notation  is 
similar  and  we  follow  closely  their  development  through  their  Eq.  (35). 

II.  TOE  INTERMEDIATE-AVERAGE  MUTUAL  COHERENCE  FUNCTION 

For  a  stationary  scalar  wave  field,  the  mutual  coherence  function 
for  the  correlation  of  two  space-time  points  is  often  written 

T 

r(P1,P2,x,T)  =2^  |  VT(P1,t+x)  V*(P2,t)  dt,  (1) 

-T 

where  the  limits  for  the  time  integration  are  allowed  to  approach  infi¬ 
nity.  For  this  case,  however,  we  wish  to  keep  the  parameter  T  finite. 


and  by  the  subscripts  indicate  that  we  assume  a  knowledge  of  V^CP^t+t) 
and  VT(»'2,t)  only  over  the  finite  sample  length  2T.  We  wish  to  call 
"Pl»P2»T»*0  *n  te  1186 di at e - ave rage  mutual  coherence  function  and 

carefully  stress  that,  for  arbitrary  T  or  shift  of  origin,  it  may  boar 
little  resent) lance  to  the  mutual  coherence  fraction  defined  by  the 
enseri>le  average - 

Following  Ref.  12,  we  represent  V(P,t)  as  the  temporal  Fourier 
transform  of  the  complex  analytic  signal 

00 

VT(P,t)  =  |  vT(P,a>)  exp(-iut)  du,  for  0  <  |t|  <  T  (2a) 

0 

and 

T 

vT(P,u>)  =  I  VT(P,t)  exp Citut)  dt.  (2b) 

-T 


S instituting  Eq.  (2a)  into  Eq.  (1),  interchanging  the  order  of  integra¬ 
tion  and  time  averaging  and  performing  the  time  integration,  we  get 


OD 

I’CJ'j.P^-r.T)  =jj  WT(P1,P2,cd1,(1)2)  exp^iujT)  sinc[ ((dj-w^T]  dt^d^,  (3) 


where 


sine  x  = 


and  the  function 


WT^PrP2*“l*“2^  ~  vt(Pi»“i)  vt^P25“2^; 


C4) 


(5) 
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and  the  subscript  T  here  and  later  implies  a  function  based  on  the  elec¬ 
tric  field  statistics  only  for  the  particular  sample  2T  in  length  about 
the  origin.  The  sine  function  of  Eq.  (3)  assumes  the  role  of  a  low-pass 
filter.  It  T  is  very  small,  the  two  frequency  variables  of  Eq.  (3)  are 
essentially  independent  and  all  cross  terms  are  represented  in  the  pro¬ 
duct  of  Eq.  (S) .  These  cross  terms  form  a  high-frequency  spectral  con¬ 
tent.  However,  as  T  tends  to  infinity,  the  sine  function  assumes  the 
role  of  a  delta  function,  constraining  correlation  to  occur  only  between 
identical  frequencies  in  the  transform  product  and  forcing  the  integral 
to  a  cne-dimensional  form.  In  the  limit  jf  large  T,  the  filtered  spec¬ 
trum  of  Eq.  (3)  becomes  the  mean  square  value  (dc)  of  each  temporal 


frequency  component  in  the  signal. 

Following  Mardtand  and  Wolf*2  and  the  earlier  lead  of  Walther,13 
vT(P,w)  is  represented  in  the  form  of  an  angular  (spatial)  spectrum  of 
plane  waves  in  Cartesian  coordinates  where 

00 

vT(P,<o)  =  ||  apfp.q.w)  exp[ik (px+qy+mz)  dpdq,  (6) 


,,  2  2.1/2 

m  =  (1  -  p  -  q  ) 


2  2 

if  p  +  q  <  1 


-  ,  2  A  2  ,.1/2 

=  x(p  ♦  q  -  1) 


2  .  2  . 
if  p  +  q  >  i. 


k  =  u/c. 


where  c  is  the  vacuum  velocity  of  light.  Equation  (6)  indicates  that 
vT(P,w)  is  formed  by  a  superposition  of  homogeneous  spatial  waves  propa¬ 
gating  in  the  half  space  z  >  0  for  the  criterion  expressed  by  Eq.  (7a) 
and  a  set  of  evanescent  waves  propagating  parallel  to  the  plane  z  =  0 
for  the  case  described  by  Eq.  (7b). 
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Expressing  Eq.  (5)  in  the  fora  of  the  angular  spectnm  of  plane 
waves  defined  by  Eq.  (6)  we  get 


WT*Pl,P2*“l*“2*  =  IIW 

— OO 

x  expfikjCPjXj  +  q^  +  m^)]  (9) 

/ 

X  e^)[-ik2(p2x2  ♦  q2y2  ♦  «2Z2)]  dp^dp^ , 

where 

AT(P1»q1;P2»‘l2;“l,<,,2)  5  aTtp2,q2;u2)*  (10) 


Using  Eq.  (9),  therefore,  the  intermediate -average  mutual  coherence 
function  of  Eq.  (3)  can  be  written 

oo 

r(Pl»P2»T»T)  =  ||  exp(-iu)jT)  sincfuj  -  u2)T]  du1dUi)2 
0 
OO 

x  Wll  »q1;p2»q2;‘01»“2)  (^) 

—00 

X  e^jfikjCpjXj  +  q^  +  Bjzp] 
x  expt-ik^^  +  w2  *  OjZ 2)]  ^dq^dq^ 


If  T  is  allowed  to  approadi  infinity,  the  limiting  form  of  the  sine 
function  forces  and  the  intermediate-average  mutual  coherence 

function  clearly  reduces  to  the  form  of  Eq.  (13)  of  Ref.  12  by  the 
elimination  of  one  of  the  time -frequency  integrals. 
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III.  THE  ANGULAR  CORRELATION  FUNCTION  AND  OTHER  SPATIAL  CORRELATION 
FUNCTIONS 

The  cross -spectral  density  function  *t(pi>P2’<‘>1,<1>2^  is  now  expressed 
as  a  four-dimensional  spatial  Fourier  integral  as 


w 


x  exp[i(flXl  +  g^  +  f2x2  +  g2y2)]  dfjdgjd^dg^ 


Equation  (9)  with  Eq.  (12)  therefore  implies 


2  2 

AT^>l,qi’P2,q2;t°l,Il52^  =  1^2  N^klpl*kiql,Zl*~k2P2,-k2q2’Z2*“l*“2^ 

x  e^I-iOCjiBjZj  -  l^m^)]  , 


(12) 


(13) 


A 


and  specifically  if  ~  z2  =  0,  .Eq.  (13)  becomes 


2  2  ~ 


AT(P1»q1:P2»‘l2;“l’“2)  =  kl  k2  WT(k1P1,klql;0;~k2P2,'k2q2;0;“l,w2)*  (14) 


Equation  (14)  indicates  that  the  angular  correlation  function 
AT^>l,ql’p2,q2’<1)l,w2^  an<^  f*16  four-dimensional  spatial  Fourier  transform 
of  the  cross -spectral  density  function  are  related  at  the  plane  z  =  0  if 

fl  =  klpl’  gl  =  klV  f2  =  _k2p2  *  OTd  «2  =  "k2q2* 

Further,  the  spatial  transform  of  the  time  spectrum  of  the  field 
can  be  represented  by 

OO 

vT(x,y,z;u)  =  jj  .vT(f,g;z;u)  exp(i(fx+gy)]  dfdg.  (15) 
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. .  «Vj»  Jkvuv  it  ttidd  t  'I  Hj  lS'.  I .  Wf!,  J  rn:  M 


Comparison  of  Eqs,  (15)  and  (6)  indicates 


^(P^w)  =  *  vT(kp,kq,z,u)  exp(-ikaz) 


Using  Eq.  (16)  in  Eq.  (1))  gives 


Aj"(Pj*<lj»P2»<l2*<,,l,w2^  =  ^1^2  ^‘r^klPl,kl^l,Zl’<l,l^ 

x  v*(k2p2,k2q2,z2;u2) 


x  ejqjf-iOCjiDjZj  -  k2m2z2)  ] . 


The  inter  ediate-average  spatial  correlation  function  is  defined 


^(klPrkiqi’2rUl)  ^T^C2P2’k2<l2,Z2;t02-) 


=  VT(k1p1,k1q1;z1;k2p2,k2q2;z2;u1>a)2) 


=  VT(k1p1,k1q1;0;k2p2,k2q2;02;u1,to2) 


>:  e3qp[i(k1m1z1  -  k^z^]  , 


(18a) 


(18b) 


where  the  product  form  of  Eq.  (18b)  is  inqilied  by  the  independence  of 
tine  left  side  of  Eq.  (17)  on  z^  and  z2.  If  we  set  z^  =  z2  *  0,  Eqs.  (17) 
and  (18)  inply 

Ar(Pl»qi»P2»<l2*,,>l,a>p  =  ki^k2  jP  ’k2p2  ,k2^2*^’^^  *^2^  * 

Finally  using  Eqs.  (14)  and  (19),  we  find 

^T^l*^l*®’^2*®2’^*ul,&>2^  =  ’  (20) 
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Ihus  the  relationships  between  the  angular  correlation  function  and  the 
cross-spectral  density  function  are  established  by  Eq.  (14)  and  the 
angular  correlation  function  and  the  spatial  frequency  correlation 
function  by  £q.  (19)  for  the  case  of  intermediate  time  averaging. 

IV.  THE  INTERMEDIATE-AVERAGE  CORRELATION  FUNCTION  IN  THE  FAR  FIELD 
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Now  the  form  of  the  cross -spectral  density  function  is  examined  in 
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the  far  field.  Defining  r^  =  (xr  ■>  yn  +  z  )  ,  we  seek  the  asymp¬ 
totic  forms  for  the  case  of  and  P 2  tending  to  infinity  in  the  paths 
indicated  by  the  direction  cosines 
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Rewriting  Eq.  (9)  using  the  definition  of  Eq.  (10),  we  have 


03 

WTCPl’P2;Ul,a>2)  =  ||  aT^l»‘*l;wl)  exPtiki^ixi  +  qlyl  4  mizi^ 


||  4(P2><l2;a)2J  exP["ik2(p2X2  +  q2X2  +  ®2Z2)''  ^2dq2  * 


As  •  and  k^  ■*  00 ,  the  asymptotic  form  of  the  two-dimensirnal 

integrals  is  given  by  Miyamoto  and  Wolri4  as 


•  v  y 

cos8i  cose2  Vr^  •  ;  ~F2  •  ij  ^i*"?3 

explif^rj  -  k2rp] 


3  -  , 


where 


z. 

—  =  cos 6,  *  —  =  cose.  (23) 

T1  1  r2  2 

and  use  of  Eq.  (10)  has  been  Bade.  Finally  if  Eq.  (22).  is  substituted 
into  Eq.  (3),  we  get  for  the  intermediate-average  mutual  coherence 
function 


r(P1,P2»x»T)  =  ||  exp(-iu1T)  sinc[(u1  -  w2)T]  du^do^ 
0 

e^p[i(k1r1  -  k2r2)] 


4w 

X  cose j  cose2 


rlr2 


(24) 


X1  yl  X2  y2 
Vrj  *  Ti  *  r2  '  r2  ;  * 


In  addition,  because  of  the  relationship  given  earlier  relating  the  angu¬ 
lar  correlation  function  to  the  cross -spectral  density  and  the  spatial- 
frequency  functions  (Eqs.  (14)  and  (19)),  Eq.  (22)  can  be  written  in  the 
following  forms: 


exp[i(k  r  -  k  r  )] 

WT(P1,P2;u1,iiJ2)  ^  4ir  cosO^  cos02  kjk2  - 


rlr2 


(25a) 


a  klxl  Vl 


Tv  r. 


•  o  *  - 
>  v  > 


k2*2 


k2y2 

—  ;  o  ; 


=  4ir  cos0j  cose2  kjk2 


eaqjfiOCjrj  -  k2r2)] 


rlr2 


(25b) 


X1  yl  *2  y2 

x  V_(kt  —  ,  k.  —  ;  0  ;  k_  ,  k.  -=■  ;  0  ;  <*>,,«.,)  . 

i  i  Tj  1  Tj  c  r2  z  i  z 
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Equations  (25a)  and  (25b)  can  also  be  used  in  Eq.  (3)  to  provide  alter¬ 
nate  fonas  of  Eq,  (24).  Equations  (24)  and  (25)  form  the  modified 
version  of  the  forms  given  in  Ref.  12,  Eqs.  (33)  through  (35).  Using 
these  results,  we  are  now  in  a  position  to  form  the  self-intensity 
function  in  the  far  field. 

V.  THE  SELF-INTENSITY  FUNCTION  IN  THE  FAR  FIELD 

We  now  examine  the  form  of  the  self-intensity  in  the  far  field  by 
letting  points  Pj  =  =  P  and  then  letting  the  time  delay,  t,  be  zero. 

Under  these  conditions,  the  mutual  coherence  function  reduces  to  the 
self- intensity, 15  and,  using  Eq.  (25b)  in  Eq.  (3)  and  the  definition 
given  in  Eq.  (18a),  we  have 


I  (P,T)  =  4ir2  cos2e  ||  dwj 


u 


dii) 


2  2 
c 


2^2  exp[i(k  -  k  )r) 

- 2 - sinc[  (tox 


u2)T] 


(26) 


*  vT(kx  f>  r  J  0  JWj)  v*(k2  k2  r  ;  0  ’“2^  * 


where,  as  indicated  earlier,  the  sine  function  acts  to  suppress  temporal 
frequencies  in  the  cross  spectrum  higher  than  ‘v  1/(2T)  Hz.  We  now  uti¬ 
lize  the  linear  transformation  of  the  time  frequency  variables  (for 
which  the  Jacobian  is  mity)  defined  by 


U1  ~  u2  H  p 


and 


“l  +  “2  _ 
— - - -  0 


(27) 


Writing  the  w  variables  in  terms  of  these  center-of-mass  coordinates 
we  get 

2o  +  p  ,  2a  -  p  _ _ 

“l  ~  — 2 — ~  and  “?  ~  — 2 — ~  ’  (28) 
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which  when  substituted  into  Eq.  (26)  gives 


2 

I(P,T)  a  (^  cos20  |  da  a2  |  sinc(pT)  exp(ipr/c)  dp 

0=0  0=-eo 


(29) 


xv;(^,^;0;a-p/2), 


where  the  dependence  of  the  amplitude  on  the  difference-frequency  co¬ 
ordinate,  p,  has  been  dropped  since  for  quasi mono chromatic  radiation 
a2  »  |p2/4|. 

Now  using  the  defining  transform  relation  of  Eq.  (15),  we  write  the 
spatial  correlation  function  at  the  source  (z=0)  where 

o  ;  -P/2)  ,  2^4  ;  0  ;  o-,m 


2  lljj  VT^1  ,r'l ’  ®  ■  — P/2)  vj(^2  ’^2  *  ^  *  P“p/2) 


C30) 


exp[-i(2i^)  (?1  f  *  r,j  i)]  exp  [i (2^)  ft 2  i  .  n,  *>]  d^dpjdr,. 


Now  the  intermediate-average  spatial-correlation  function. 


vT(C1,n1;  0  ;  a+p/2)  v*(C2,n2;  0  ;  a-p/2) 


when  considered  with  the  filtering  action  of  the  sine  function  of  Eq.(29) 
will  have  an  effective  contribution  only  for  the  low  frequency  components 
formed  by  the  difference -frequency  terms  1/(2T)  .Hz  or  less. 


18 


In  addition,  we  assume  the  node  population  to  be  a  slowly-varying 
function  of  a,  since  a  »  p/2  and  thus  we  write 

ty^n^Oja+e/2)  v*(?2,n2>0;o-p/2)]low  freq> 

*  A(Cj ,n1;(m>/2)  A(52,n2;o-p/2)  expfi.KSj  -n2;p)  (31) 

-  A(j-.,o)  A(£2,o)  H(p)  expfi*^  -i2:p)], 
where 

(  1  for  A(£;c±p/2)  /  0 

H(p)  =  J 

|  0  otherwise. 

Equation  (31)  acknowledges  the  loss  of  the  optical-frequency  phase, 
while  maintaining  the  phase  of  the  intensity  envelope  formed  by  temporal 
beat  codes.  The  degree  to  which  the  phase  of  this  envelope  is  detected 
depends  on  the  bandwidth  of  the  source  and  the  detector  resolution,  2T. 

Essentially,  these  arguments  ware  made  by  H anbury  Brown  and  Tfciss1® 
except  for  the  defining  of  the  H  function.  Its  introduction  is  brought 
about  by  the  description  of  narrow-band  sources  b>  terns  xn  A(£,o).  For 
a  thermal  source  of  relatively  large  bandwidth,  the  maximum  difference 
p  will  extend  far  beyond  th<  tenoral -frequency  response  of  the  system 
(here  reflected  in  the  sine  tern  of  Eq.  (24))  and  be  continuous  as  well. 
But  for  a  lnsor  source  exhibiting  a  series  of  axial  modes,  the  complete 
diffexvnce-fit'qufrncy  duaain  might  lie  entirely  within  the  system  response 
bet  be  pi?ce-wise  continuous  in  its  extent. 

Relative  to  the  representation  of  the  intermediate-average  by  the 
form  of  Eq.  (31),  we  wish  to  reiterate  a  statement  made  following  Eq.(l) 
that  the  intermediate-averaging  process  may  bear  little  resemblance  to 
the  infinite  tiwe  average,  even  so  far  t«s  the  detail  cf  the  aim  Jitudc 


terms,  A (£,o).  This  situation  would  be  serious  if  our  intent  were  to 
infer.,  for  example,  the  time-frequency  statistics  of  the  source.  But  in. 
the  present  concept,  we  desire  only  to  infer  the  spatial  properties  of 
the  source.  If'  we  consider  a  multi-axial-mode  laser  beam  scattered  from 
a  spatially-rough  surface,  the  lack  of  correspondence  between  the  two 
average  is  uninportant,  for  all  such  mode  history  is  integrated  out;  all 
areas  of  the  scatterer  see  the  same  mode  characteristics.  Any  mode 
fluctuation  would  be  seen  simply  as  a  variation  in  total  received  power 
from  one  sasple  to  the  next.  Here,  we  simply  require  for  one  detector- 
resolution  time  over  a  spatial  domain  that  the  process  of  Eq.  (31)  ex¬ 
hibit  a  minimum  of  two  temporal  modes  (to  maintain  the  phase  term 

-  J»2»p}  “i**1  sufficient  mode  population  (reflected  in  the  amplitude 
terms  A(£,o))  such  that  quantum  noise  in  both  the  carrier  wave  and  the 
detector  can  be  ignored. 


Using  the  results  of  Eq.  (31)  in  Eq.  (29)  and  taking  0  «  1,  we 


write 


I(P»T)  -  (^)  |  2  dc  |  dp  exp(ipr/c)  sinc(pT)  H(p) 

0  -«> 

x  jjjj  -  £-  ;?)3  C32) 

*  f  *  Oj  *>]  ejp[i(£^-2-Ke2  f  «  r,2  *)]  SSjdtjdnj*,. 


We  now  have  the  self  intensity  in  the  far  field  expressed  as  a  dothle 
integral  over  sum-  and  difference -frequency  conponents  as  well  as  two, 
two-dimensional  spatial  Fourier  transforms  over  the  source.  We  will 
next  use  it  to  form  the  fourth-order  correlation  function. 


We  now  investigate  the  phase  variation  in  the  above  terns .  In  order 
to  develop  the  condition  of  spatial  incoherence,  we  demand  that  the  source 
surface  be  spatially  rough,  specifically  that  phase  changes  of  m  occur 
'b  10~6  m  across  its  extent.  Next  we  make  an  order  of  magnitude  calcula¬ 
tion  of  the  size  of  the  term  g~p^2  (— )  £  ^  _  (Soosing  some  approx- 

c  r  cr 


.  .  cj  2ir 

imate  values,  —  =  r— 
c  X 

n.  103  m"1.  Lettinj 


—  =  — ■  n.  10  ^  m"1,  x  *b  lm,  r  'v*  10*  m,  we  find 

C  A 

3  -3 

Letting  ir  =  10  £  i  up  lies  K  must  move  through  ir  x  10  m 


to  get  ir  change. 

Thus  the  phase  terms  of  Eq.  (33)  describing  the  source  spatial 
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roughness  have  a  spatial  rate  of  oscillation  of  about  10  greater  than 
the  remaining  terms.  On  the  basis  of  a  somewhat  heuristic  argument 
given  in  the  Appendix,  we  assert  that  the  dominant  terms  of  Eq.  (33)  are 
those  for  which  the  phase  variation  is  a  minimum. 

Upon  examining  Eq.  (33),  we  note  that  in  the  first  cosine  argument, 
the  random  phase  terms  add  and  thus  there  is  no  choice  of  £  or  p  terms 
which  will  slew  the  phase.  Turning  to  the  second  cosine  expression,  we 
see  that  there  are  two  general  cases  for  which  the  phase  variation  is  a 
minimum.  The  first  is  the  familiar  one  for  which  the  parameter  p  is 
equal  to  zero.  This  corresponds  to  the  normal  -.nfinite-time  average 
situation.  For  this  case  we  set  and  =  j^'  with  the 

result  that  all  the  source-spatial  information  collapses  in  the  cosine 
argument.  This  is  the  normal  dc- intensity  term  in  the  far  field.  For 
convenience,  we  will  suppress  this  term  in  following  expressions.  For 
P  0,  the  minimum  phase  can  be  achieved  by  removing  the  primes  of 
Eq.  (33),  halving  the  nunfoer  of  integrals  to  give 


Re  I(P  ;T)  Re  I(P  ;T)  =  ( 

C«>  oio 


4  , 
a  da 


(O 

j  dp 


sine  (pT)  H(p) 


X  ||  I(£2,0))  expj-i|  [y* 

oo  V 


x2)  ♦  n2(yry2)  }d52dn2 


dgjdnj 


(38a) 


']) 

||  i(ir«)  exp|-i~  [wv  +  ni(yry2)]|j 


— — j  [  H(p)  sinc^(pT)  dp  f  w*  dto 

fcrl4  J  J 

lcrj  p-0+  01=0 


||  1(5,0.)  exp  j-i|  ^(Xj-Xj)  +  n(yry2)j|d5dn 


where  the  mean  intensity  has  been  suppressed. 

If  the  source  time-frequency  characteristics  of  the  source  are  such 
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that  p  is  continuous  over  the  domain  for  which  sine  (pT)  is  essentially 
nonzero,  then  H(p)  =  1  for  all  p  and  then  the  difference-frequency  term 
can  be  integrated  easily  to  give 


(38b) 


sine  (pT)  dp  =  ir/(2T). 


(33) 


p=0+ 


It  is  evident  from  Eq.  (38b)  that  the  two-point  intensity  corre¬ 
lation  in  the  far  field  is  proportional  to  the  modulus  of  the  spatial 
Fourier  transform  across  a  spatial ly-rough  source.  Only  for  the  case 
that  the  intensity  distribution  on  the  source  is  pure  even  can  the  phase 
of  the  spatial  transform  be  inferred  and  used  to  invert  uniquely  Eq.  (38) 
to  get  the  intensity  distribution  on  the  source,  I (£ ,u) . 
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VII.  SPATIAL  AVERAGING 


Let  us  rewrite  Eq.  (38)  in  a  position-vector  notation  for  the 


soatial  coordinates  so  that 


c* 

I(PrT)  I(P2,T)  =  |  H(p)  sinc‘(pT)  j 

^•Cl^  p=0+  u=0 


f  ^ 


|  1(1,05)  exp(-i  £i-d)  dl  , 


where 


-  s  (xrV  1  +  0^i-y2)  J  » 


1  =  5i  +  nj 


By  observing  Eq.  (40)  one  can  define 


KPjjT)  I(P2;T)  =  Cdj.^jT)  =  C(d;T) 


It  is  thus  obvious  that  the  two-point  intensity  correlation  for  the 
ideal  case  of  noiseless  imaging  is  spatially  dependent  only  on  the  co¬ 
ordinate  difference  in  the  far  field.  This  property  of  Eq.  (40)  reflects 
its  intrinsic  spatial  stationarity,  indicating  that  the  global  dependence 
is  absent.  If  noise  introduced  in  the  imaging  process  is  independent  of 
the  source  character,  then  it  will  be  spatial  in  character  and  related  to 
atmospheric  turbulence  or  detector  limitations.  Averaging  the  function 
of  Eq.  (40)  over  an  area  woul'’  increase  the  signal-to-noise  ratio,  the 
area  of  sufficient  size  being  determined  by  the  correlation  interval  of 


the  noise  statistics. 


I  j 


Averaging  Eq.  (40)  over  an  area  Aq  singly,  then,  multiplies  the 
two-point  correlation  by  a  factor  Aq  giving 


<  I(Pt;T)  I(P_;T)  >  _.  .  =  <  C(d;T)  >  . 

v_i»  *  —2  3  spatial  average  —  Aq 

»  QO  00  / 

— ~  [  H(p)sinc2(pT)  dp  f  a>4  dw  f  1(5,0))  exp(-i£  £*d)  d£ 

(Crv  3  j  )  T 

'  J  0=0+  «1)=0  -CO 


VIII.  SUMMARY  AND  CONCLUSIONS 

Using  a  methodology  of  Marchand  and  Wolf  as  a  starting  point,  we 
have  developed  an  expression  for  the  two-point  intensity  correlation  in 
the  far  field,  independent  of  time  averaging  except  for  the  temporal 
resolution  of  the  detector.  This  result  is  derived  for  narrow-band,  high- 
intensity  light  scattered  from  a  spatially-rough  surface. 

To  illustrate  the  above  ideas,  we  wish  to  describe  a  simple  experi¬ 
ment  which  embodies  these  mathematical  ideas.  Using  a  helium-neon  laser 
in  a  single-axial-mode  configuration,  a  symmetrical  source  is  transil- 
luminated.  The  source  has  a  random-phase  character  to  obtain  spatial 
incoherence.  The  far- zone  intensity  pattern  is  then  recorded  by  film 

using  an  exposure  time  less  than  the  reciprocal  of  the  source 
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bandwidth.  This  insures  that  the  phase  of  the  beat  frequency  is 
recorded. 

Next,  the  film  is  developed  so  that  it  is  linear  in  intensity  and 
used  to  make  two  identical  positive  transparencies.  The  positives  are 
then  placed  in  a  collimated  beam,  forming  the  correlated  intensity  over 
an  averaging  area,  Aq.  The  signal  transmitted  by  the  transparency  pair 
is  optically  Fourier  transformed,  a  dc  block  inserted  to  remove  the  ni- 
wanted  average  term,  and  the  total  remaining  irradiance  measured.  This 
signal  represents  the  mathematical  expression  given  by  Eq.  (44)  for  the 
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transparency  spatial  lag,  d.  Since  the  source  is  known,  a  priori,  to  be 
symmetrical,  the  transform  of  Eq.  (44)  is  pure  real.  The  square  root  of 
the  correlation  signal  is  proportional  to  the  spectrum  which  is  then 
known  as  a  function  of  spatial  lag  <1.  Finally,  this  one  dimensional 
signal  is  Fourier  transformed  by  machine  to  give  the  scaled  source 
irradiance. 

We  have  therefore  shown  that,  given  a  symmetrical,  spatially 
incoherent  source  illuminated  by  high-intensity  light,  the  far- zone 
intensity  pattern  can  be  used  to  form  the  optical  image  of  the  source  if 
the  signal  is  recorded  with  suffj  ciently-short  time  resolution. 


APPENDIX 

We  argue  heuristically  that  the  form  of  Eq.  (33)  for  arbitrary 
<KjL,  “  Ka  IP)  approaches  a  condition  expressed  by  one  form  of  the 

^  JO 

Riemann-Lebesgue  theorem  which  asserts  that 


U 

lim  F(x)  cos(kx)  dx  =  0  . 
k  -*•  »  J 


(Al) 


Here,  tlie  terms  on  the  right-hand  side  of  Eq.  (33),  independent  of 

♦(5j  -  JS^JP)*  assume  the  role  of  F(x),  while  cos(kx)  is  represented  by 

terms  dependent  on  +(Cj  -  j^JP)  *  Taking  the  limit  in  Eq.  (Al)  in  this 

context  implies  a  surface  roughness  of  arbitrarily  small  scale  and  that 

the  energy  will  be  radiated  at  these  very  small  spatial  frequencies . 

These  waves,  though*  are  evanescent  and  do  not  transfer  any  energy. 
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This  kind  of  contradiction  is  discussed  by  Be  ran  and  Parrent  and  can 
arise  similarly  when  the  constraint  of  extreme  spatial  incoherence  is 
expressed  in  the  form  of  a  delta  function  of  coordinate  differences  and 
interpreted  in  a  strict  sense.  On  the  basis  of  the  limit  implied  by 
Eq.  (Al),  we  infer  that  the  dominant  terms  of  Eq.  (33)  are  those  for 
which  the  phase  variation  is  a  minimum. 


This  similarity  to  second-order  spatial  incoherence  can  be  realized 
by  a  physical  argument.  If  we  consider  the  imaging  of  a  spatially- rough 
surface  illuminate  *  *»y  perfectly  monochromatic  light  under  the  condition 
that  the  surface  i  sufficiently  large  to  be  composed  of  many  independent 
scatterers,  we  ■  ~  see  that  the  light  waves  from  each  scatterer  will  be 
uncorrelated  at  a  point  at  the  receiver  by  virtue  of  the  random  phase 
introduced  at  the  source.  Further,  the  intensity  at  this  detector  point 
will  simply  be  the  sum  of  the  squares  of  the  individual  field  contribu¬ 
tions  at  that  point.  It  should  then  seem  reasonable  that  if  we  set  tip 
a  similar  situation  of  random-phase  addition  in  difference-frequency 
space,  the  incoherent  addition  of  intensities  should  hold  as  well,  without 
appealing  to  the  approach  of  extending  an  average  over  a  time  history  of 
temporal  fluctuations. 
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